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Description 



METHOD OF SELECTING A CRITICAL 
PLANE FOR MULTI-EVENT FATIGUE LIFE 
PREDICTION 

Background of Invention 
[0001] i. Field of the Invention 

[0002] The present invention relates to a method of selecting a 
critical plane for use in predicting fatigue life of an object 
experiencing multiple stress inducing events. 

[0003] 2. Background Art 

[0004] Uniaxial fatigue theory is a fatigue life prediction tool. It is 
used to predict the ability of an object to withstand re- 
peated exposure to stresses before cracking. The process 
is referred to as "uniaxial" because the prediction is based 
on one-dimensional stresses being induced within the 
object. 

[0005] Critical plane theory is a method to model multi- 
dimensional stresses as one-dimensional stresses. The 



one-dimensional stresses can then be analyzed using uni- 
axial fatigue theory to predict fatigue life. The one- 
dimensional stresses are determined by selecting a critical 
plane direction and determining how much stress the 
multi-dimensional stress exert in the direction of the crit- 
ical plane. 

[0006] a problem with the critical plane approach arises when 
the object experiences a number of stress inducing 
events. In the past, if there were three events under con- 
sideration, three damage calculations would occur relative 
to three different critical planes. The fatigue life would 
then be determined by adding together the three damage 
calculations. This is a problem because the damage calcu- 
lations are relative to three different critical planes. Ac- 
cordingly, there exists a need for a method of selecting a 
critical plane so that the critical plane can be used for the 
damage calculations of all events. In this manner, damage 
calculations from each event can be added together for 
the purposes of calculation fatigue life using uniaxial fa- 
tigue theory. 
Summary of Invention 

[0007] The present invention overcomes the foregoing problem 
with a method for selecting a critical plane. The method 



analyzes stress induced by multiple events to select a crit- 
ical plane. The critical plane is suitable for use in calculat- 
ing damages caused by each event. The present invention 
thus eliminates the problem of adding together damage 
calculations made relative to different critical planes. 
[0008] one aspect of the present invention relates to a method 
for use with an object experiencing multiple stress induc- 
ing events. The method comprises determining a stress 
distribution histogram for each event and superimposing 
each stress distribution histogram to create a superim- 
posed stress distribution histogram. A critical plane is de- 
termined based on the superimposed stress distribution 
histogram. 

[0009] one aspect of the present invention relates to a system for 
predicting fatigue life of an object experiencing multiple 
stress inducing events. The system comprises a computer 
which determines a stress distribution histogram for each 
event. The stress distribution histograms are superim- 
posed to create a superimposed stress distribution his- 
togram so that the computer can select a critical plane 
from the superimposed stress distribution histogram. 

[0010] one aspect of the present invention relates to a method 

for predicting a fatigue life of an object experiencing mul- 



tiple stress inducing events. The method comprises deter- 
mining a stress distribution histogram for each event and 
weighting at least one of the stress distribution his- 
tograms. The method further comprises superimposing 
each stress distribution histogram to create a superim- 
posed stress distribution histogram. A critical plane is cal- 
culated based on the superimposed stress distribution 
histogram. Once that critical plane is calculated, the 
stresses induced by each event are adjusted relative to the 
critical plane. The largest of the adjusted stresses is de- 
termined. New events are determined by shifting each 
event according to the largest stress. This provides mean 
stress correction. The damage caused by each new event 
is calculated. The damage associated with the largest 
stress is subtracted from each of the damage calculations 
to determine a corrected damage for each of the new 
events. The corrected damages are added together to de- 
termine a total damage. The total damage is used to pre- 
dict the fatigue life of the object. 
Brief Description of Drawings 

[° 011 ] Figure 1 illustrates a system in accordance with the 

present invention; 
[0012] Figure 2 illustrates a two-dimensional stress sensor; 



[0013] Figure 3 illustrates a stress distribution histogram for an 
event #1; 

[0014] Figure 4 illustrates a stress distribution histogram for an 
event #2; 

[0015] Figure 5 illustrates a stress distribution histogram for an 
event #3; 

[0016] Figure 6 illustrates a superimposed stress distribution 
histogram; 

[0017] Figure 7 illustrates a weighted stress distribution his- 
togram for the event #1; 

[0018] Figure 8 illustrates a weighted stress distribution his- 
togram for the event #2; 

[0019] Figure 9 illustrates a weighted stress distribution his- 
togram for the event #3; 

[0020] Figure 10 illustrates a superimposed stress distribution 
histogram for the weighted stress distributions his- 
tograms; 

[0021] Figure 11 illustrates a method for predicting a fatigue life 
of the an object experiencing multiple stress inducing 
events; 

[0022] Figure 12 illustrates a stress history for multiple stress in- 
ducing events; 
[0023] Figure 13 illustrates stress history of one event; 



[0024] Figure 14 illustrates a portion of the stress history shown 
in Figure 15 corresponding with a largest stress; and 

[0025] Figure 15 illustrates a new event determined in accor- 
dance with mean-stress correction. 
Detailed Description 

[0026] Figure 1 illustrates a system 10 in accordance with one of 
the many aspects of the present invention. The system 
generally comprises a computer 16, an object 18, and a 
transducer 20. The object 18 experiences multiple stress 
inducing events. The computer 16 operates in conjunction 
with the description continued below to predict the fa- 
tigue life of the object 18. The prediction is based on a 
critical plane approach. The critical plane approach of the 
present invention selects a single critical plane for use in 
calculating damages caused by each event. Accordingly, 
the damages caused by each event are calculated with re- 
spect to the same critical plane and can be added together 
to predict fatigue life. 

[0027] while not intended to limit the present invention, the ex- 
amples described herein relate to an automotive testing 
environment. The automotive testing environment sub- 
jects automotive vehicles to multiple stress inducing 
events. Each event induces different multi-dimensional 



stresses. The present invention can predict the fatigue life 
of the object with a minimal amount of data regarding the 
multi-dimensional stresses induced by each event. The 
prediction can be used for any number of computations, 
including to design a vehicle component according to de- 
sirable performance criteria. 

[0028] Automotive testing generally comprises driving a vehicle 
over a number of simulated road conditions. The road 
conditions simulate driving conditions the vehicle may ex- 
perience in the future. The performance of the object 18 
under the test conditions can provide valuable feedback 
as to how well the object 18 has been designed to with- 
stand future operation. 

[0029] Each different road test induces stresses differently within 
the object 18. Accordingly, each road condition is consid- 
ered as a separate multi-dimensional stress inducing 
event. The transducer 20 measures deformation induced 
by the stresses in the object 18. The computer 16 compu- 
tates the multi-dimensional stresses induced by the stress 
inducing event from the deformation. The effect of the 
stresses on fatigue life is calculated from one run through 
each event rather than cycling through each event until a 
crack develops. This approach saves time, computational 



costs, and allows components to be modeled with less 
testing time such that further re-designs can be com- 
pleted with less design revisions. 

[0030] Figure 2 illustrates the stresses calculated by the com- 
puter 16. As shown, normal 24 and shear 26 stresses are 
calculated. The calculated stresses are two-dimensional. 
Only these stresses are calculated because the present in- 
vention makes a number of assumption that reduce typi- 
cal complexity of such calculations. Namely, the typically 
used nine component stress tensor is reduced to the nor- 
mal stress 24 and shear stress 26. The reduction is done 
by assuming (i) the object 18 is relatively thin in compari- 
son to its longitudinal and latitudinal dimensions, (ii) that 
the stresses occurring throughout the thickness are ig- 
norable as they are relatively insignificant compared to 
the other stresses, and (iii) that the stresses occurring on 
opposite sides are equal. However, other assumptions can 
be made to include more or less stress calculations. 

[0031] Figures 3-5 illustrate histograms for three events referred 
to as event #1, event #2, and event #3, which are respec- 
tively referred to with reference numerals 34, 36, and 38. 
The histograms are a means to track the stresses induced 
by each event. The histograms are generated by the com- 



puter. The computer breaks up each event into a number 
of time intervals. The stresses (normal stress 24 and shear 
stress 26) induced by each event are calculated for the 
time period. A principle stress and an angle of the princi- 
ple is used to represent the stresses as a one-dimensional 
stress. The principle stress represents in one-dimension 
the cumulative effect of the normal stress 24 and the 
shear stress 26. The angle represents a direction of the 
principle stress. 

[0032] The histograms utilize a binning system to track occur- 
rences of the principle stresses and angles over the period 
of time that the object experiences stresses. The principle 
stress must have an average value within a bin for a pre- 
determined period of time to count as an occurrence. The 
binning system organizes the principle stresses along a 
vertical axis and the angles along a horizontal axis. The 
vertical axis includes a range from 0 to 100 MPa and 0 to 
-100 MPA. The ranges span positive and negative bins of 
0-10, 10-20, 20-30, 40-50, 50-60, 60-70, 70-80, 
80-90, and 90-100 MPA. The horizontal axis includes a 
range from 0-90° and 0 to -90°. The ranges span positive 
and negative bins of 0-10, 10-20, 20-30, 30-40, 50-60, 
60-70, 70-80, and 80-90. Please note, the angles and 



principle stresses are relative to the center of the stress 
plane shown in Figure 2. 

[0033] Figure 6 illustrates the superimposing of the histograms 
shown in Figures 3-5 to create a superimposed histogram 
40. The single critical plane is selected from the superim- 
posed histogram. This overcomes the problem of adding 
together damages calculated relative to different critical 
planes by determining a single critical plane for all the 
events. The damages caused by each event can then be 
determined with respect to the same critical plane. Dam- 
ages calculated with respect to the same critical plane can 
be added together to predict fatigue life. 

[0034] The superimposed histogram 40 is formed by adding to- 
gether the occurrence from each event that fall within 
common bins. For example, in exemplary bin 42, event #1 
includes 63 occurrences of a principle stress with an am- 
plitude in the range of 0-10 MPA and an angle in the 
range of 40-50 degrees, bin 42' of event #2 similarly in- 
cludes 35 occurrences, and event bin 42 M of #3 similarly 
includes 0 occurrences for a superimposed bin 42 m total 
of 98 occurrences. The same process is repeated for each 
bin to populate each bin in the superimposed histogram 
40. The single critical plane is then selected from the su- 



perimposed histogram 40 by a computer or other device. 
This allows the effects from each event to be included 
when selecting a single critical plane. 

[0035] Figures 7-10 illustrate an aspect of the present invention 
in which the histograms are weighted prior to superim- 
posing the histogram. The weighting is done so that some 
events can affect fatigue life more than the other events. 
The weighted histograms of events #1 to #3 are respec- 
tively referred to with reference numerals 46, 48, and 50. 
For exemplary purpose, it is assumed that event #1 and 
event #3 are run five times for each time event #2 is run. 
Rather than run each event #1 and #3 over and over 
again, the present invention multiples a single run of the 
event with an event multiplier 54. 

[0036] As shown in Figure 7, the event multiplier 54 of event #1 
is five as event #1 is to be run five times. As shown in Fig- 
ure 8, the event multiplier 54' of event #2 is one as event 
#2 is to be run only once. As shown in Figure 9, the event 
multiplier 54" of event #3 is five as event #3 is to be run 
five times. The event multiplier 54 can be the same or dif- 
ferent for each event. The event multiplier 54 weights the 
importance of each event to correspond with the need to 
simulate some events more than others. 



[0037] Figure 10 illustrates weighted superimposed histogram 
56. The weighted superimposed histogram 56 is deter- 
mined by adding the occurrences of each principle stress 
magnitude after application of the event multiplier 54. 
The critical plane is then calculated based on the weighted 
superimposed histogram 56. 

[0038] Figure 11 illustrates a method for predicting a fatigue life 
of an object experiencing multiple stress events in accor- 
dance with one aspect of the present invention. The 
method comprises a step 60 to determine a stress distri- 
bution histogram for each event. A step 62 superimposes 
each stress distribution histogram to create a superim- 
posed stress distribution histogram. A step 64 calculates 
a critical plane based on the superimposed stress distri- 
bution histogram. 

[0039] a step 66 adjusts the principle stresses shown in the his- 
tograms for each event relative to the direction of the crit- 
ical plane calculated in step 64. Once the principle 
stresses for each event are computed relative to the same 
critical plane, the damages caused by each event can be 
calculated and added together to determine the fatigue 
life of the object. 

[0040] p r j or to calculating the damages, and i n accordance with 



another aspect of the present invention, the precision of 
the life expectancy calculation is improved at a step 68 
with mean stress correction. The mean stress correction is 
typically done in accordance with the Smith-Wat- 
son-Topper rule, as understood by one of ordinary skill in 
the art. Mean stress correction is shown with more detail 
in Figures 12-15. 

[0041] Figure 12 illustrates a history graph 70. The history graph 
70 illustrates all the principle stresses for each of the 
events after adjustment to the critical plane direction. The 
vertical axis is MPa and the horizontal axis is a cycle ref- 
erence number. In the real world the vehicle is already 
prestressed, therefore, the prestressing produces a non- 
zero mean stress environment. As such, each induced 
stress has a non-zero mean stress. This can be a prob- 
lem. The problem occurs because fatigue life calculations 
determined under non-zero mean stress tend to be in- 
compatible with fatigue life calculations based on zero 
mean stress. This problem is ameliorated with mean 
stress correction. 

[0042] Mean stress correction comprises identifying each event 
within the history and shifting each event according to a 
largest stress portion 72 found the history. Once shifted, 



the events become new events. Damages are then calcu- 
lated based on the new events. 

[0043] The mean stress correction is understood by one skilled in 
the art and can be generally understood with respect to 
event #1. Figure 13 illustrates a portion 74 of history 
graph 70 that corresponds with event #1. Figure 14 illus- 
trates the largest stress portion 72. Mean stress correc- 
tion of event #1 is shown in Figure 15 and occurs when 
the event #1 portion 74 is shifted according to the largest 
stress portion 72. A new event #1, which is referred to 
with reference numeral 76, is created. The same shifting 
occurs for each event. The damages are then calculated 
from the new events such that mean stress correction can 
be incorporated and used to improve the precision of the 
damage calculations. 

[0044] Returning to Figure 13, a step 80 calculates damages 
caused by each new event according to the Palmgren- 
Miner Rule. Because the stresses were adjusted relative to 
the same critical principle plane, the damage calculations 
are in the same critical plane direction. Accordingly, the 
damage caused by each event can be added without the 
problems experienced by adding damages having differ- 
ent critical planes. 



[0045] The damages caused by each new event are inflated from 
adding the largest stress to each new event. Accordingly, 
the portion of the damages attributable to the largest 
stress portion 72 must be eliminated. The method in- 
cludes a step 82 for calculating the damage caused by the 
largest stress portion 72. A step 84 is also included for 
subtracting the damage calculated for the largest stress 
from the damages calculated for each of the new events to 
determine corrected damage for each of the new events. 

[0046] The method includes a step 80 for calculating the total 
damage from each new event based on the corrected 
damages. The total damages are calculated by adding to- 
gether each of the corrected damage calculations. This 
can be done since the damages are calculated in the same 
critical plane direction and with the same mean stress 
correction. From the total damage, a step 88 calculates 
the fatigue life. The life is calculated according to Palm- 
gren-Miner Rules by taking the inverse of the total dam- 
ages. 

[0047] The computer 16 can be used to implement the method 
according to Figure 13. More specifically, the computer 
can include a computer-implemented software program 
for predicting a life expectancy of an objected subjected 



to multiple stress inducing events as described above. The 
software program can be stored on a disc or other 
medium and loaded to a computer or other device for use 
in predicting life expectancy. 
[0048] while the best mode for carrying out the invention has 
been described in detail, those familiar with the art to 
which this invention relates will recognize various alterna- 
tive designs and embodiments for practicing the invention 
as defined by the following claims. 



